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This tutorial article treats certain aspects of the description and measurement of 
oscillator stability. Topics covered are time and frequency deviations , Allan variance , the 
zero-crossing counter measurement technique , frequency drift removal , and the “ three- 
cornered hat. ” 


I. Deviations in Phase, Time, and Frequency 

The purpose of this article is to define the Allan variance, 
relate it to various kinds of disturbances in oscillators and 
measurement systems, and to describe, from the point of 
view of an analysis and software person, the frequency stabil- 
ity measurement method used at the JPL Frequency Standards 
Laboratory. This is an expanded version of a talk given in Oct 
1983 to DSN operations. It outlines the JPL frequency stabil- 
ity measurement methodology. 


In other words, t t is the time at which the output of a perfect 
oscillator, sin (2 7tv Q t), would have the same total phase as 
oscillator i at time t. 

The relative phase and time deviations of the pair of oscil- 
lators are defined by 

m = -0 2 (o 


Let oscillators 1 and 2, both running at nominal frequency 
p 0 , have the outputs 

sin (2nv Q t + 0, (f)), sin (2 my + 0 2 (O) 

where 0,(0 and 0 2 (O are the phase deviations relative to fre- 
quency p 0 . Suppose that both oscillators are used for driving 
clocks. At true time t, the clock times r, , t 2 are defined by 


( 1 ) 

( 2 ) 


*(0 = f, - t 2 

Subtracting (2) from (1) gives 


- % 


The fractional frequency deviation is defined by 



dtp 

~di 
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2 my, = 27riy + 0 1 (O 

2 my 2 = 2 "V + 0 2 (O 


Then 


and the RMS deviations 


y(t) 


dx 

dt 




x(t) = / y(t) dt 


The average fractional frequency deviation over the time 
interval (t-r, t ) is defined by 
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where A r is the backward first difference operator. We call r 
the averaging time. When someone says that two oscillators 
differ in frequency by 5 X 1 0~ 1 4 , say, he is most likely quot- 
ing a value of y(t,r ). To understand the quote, the value of 
r must be known. Is it one second or one day? 

Figure 1 shows plots of x(?) and y(tj) from a test of two 
hydrogen masers, where r = 4094 s and both* and y are sam- 
pled once every r. The lower curve y is approximately the 
derivative of the upper curve x. Notice that x = 0 at both end- 
points. This is an artifact of the presentation; the mean y has 
been subtracted off. 

II. Allan Variance 

The Allan variance was invented to solve the problem of 
characterizing the RMS frequency deviations of a pair of oscil- 
lators. Given time-deviation samples jt(0), jc(r), •••,; c(mr) over 
the time T = mr, let 

y j = y(h, r) 

the average frequency between (/ - 1 )r and /r. One could 
compute the averages 

i m 

<*> = — “t Y x (/r) 
m + 1 ^ v 7 
/=o 

m 

iy) = — y y. 
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Often, however, these quantities tend to grow with T\ more- 
over, they may fail to be reproducible from one T-interval to 
the next. (For example, if y is a random walk, then the mean 
and standard deviation of y^ MS are both proportional to T.) 

For many types of oscillator phase deviations, the Allan 
variance, denoted by o^(r), avoids these problems while still 
yielding a meaningful measure of frequency stability. To 
define theoretical Allan variance from the preceding setup, we 
must assume that x(j't) is available for all non-negative inte- 
gers/. Then, by definition, 

1 m 

a y (T) = Hm 2(m~^T) £ 0' -T / _ 1 ) 2 (3) 

v i=2 

if the limit exists. In other words, o y (r) is the RMS average of 
the quantities ( y / - - y^ x )/y/2; thus o y (r) measures the RMS 
change iny(f,r) when t changes by r. Note that 

y f - y h j = 7 [ a 7(/>) - \x({j - i ) r)] 

= j r A 2 T x(}T) 

= i[ JC (/r) - 2 x((j - l)r) + x({j - 2)t)] 

where Aj is the backward second difference operator. The 
quantity a^(r) is called the Allan variance for the averaging 
time r. It seems that o y (r) (without the square) should be 
called Allan deviation. In this writeup, it will be done so for 
the sake of precision. In casual discussions, however, o y {r) is 
called the Allan variance without confusing anyone. 

The factor 2 (or \/2) is there for historical reasons. On the 
other hand, if the two oscillators are judged to be of like qual- 
ity, then the Allan deviation of the pair is sometimes divided 
by \Jl to give the Allan deviation of the individual oscillators. 
This is a different y/2. Moreover, this procedure assumes that 
the second r-differences of phase of the two oscillators are 
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orthogonal random processes. If linear frequency drift is 
present, or if the two oscillators are subject to the same envi- 
ronmental fluctuations, then this assumption is unrealistic, 
at least for large r. 

The usual estimator of o 2 (t), given x(f) for 0 < / < T - mr , 
is 


S 2 y (T,m) 


1 


2 (m - 1) 


E - -W 2 


/= 2 


1 

2r 7 (m - 1) 


m 2 

E [ A t*U T >] 

/= 2 


(4) 


Figure 2 shows 5 y (r, m) (“sigma”) versus r for the same 
test that yielded Fig. 1. 

A. Naive Error Estimate 

If, in the last sum of (4), the A 2 t x(jt) were independent, 
zero-mean Gaussian random variables with the same variance, 
then S 2 (r,m) would be proportional to a x 2 variable with 
m - 1 degrees of freedom. The assumptions of independence 
and zero mean are almost never realistic. Not knowing the 
true situation, however, one usually estimates the standard 
deviation of S 2 (t 9 m) by the naive formula 5^(r, m)e, where 
e = y/2/(m - 1). Then 

|y t, m) Vl - e, S y (j, m) \/l +e J (5) 


Table 1 gives a numerical example, a fragment from an actual 
test of two hydrogen masers. Elapsed time t is given in units of 
r 0 = 256 s, and x(r) is given in units of 10“ 14 s. From these 
data, the following estimates can be calculated: 


is presented as the roughest sort of “one-sigma” err or bar for 
S y {r , m). If drift is removed (Section IV), then e = \l2\(m - 2) 
is used. Figures 2 and 6 show the intervals (5) and also give 
m - 1, the “number of samples.” 
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= 8.37 X 10 


,— 16 


There are other estimators of o y {r). The RMS average of 
the second differences of x can include values of A *x(t) for t 
not a multiple of r. In the last example, for r = 3r 0 , the values 
A \x{t), t = 6, 7, 8, can be used for computing the estimate 

*5 — I4(91 2 + 87 2 + 63 2 )1 1/2 = 7.48 X 10~ 16 
n/2 3t 0 L3 v J 

Is this a better estimate than -Sy (3r 0 ,2)? If so, then in what 
sense is it better? This subject is controversial (Refs. 5, 6, 
and 8). 1 


Don’t use (5) unless m > 4; even then, don’t take it seri- 
ously. Naturally, if you know the properties of the process 
x(t) in advance, you can compute more accurate variance esti- 
mators and confidence intervals (Refs. 2, 3, 5, 6, 7, and 8). 1 
But if you did know all this, you wouldn’t be testing the 
oscillators in the first place. In any case, make sure that the 
user of the results knows the number of samples. 

B. Deterministic Examples 

The phase fluctuations of oscillators can often be modelled 
as a simple nonrandom function of time, plus a random com- 
ponent. Therefore, it is useful to know the effect of certain 
deterministic phase functions on Allan variance. 

(i) Constant phase and frequency offsets : 

x(t) = a 0 + flj t 

Here, A 2 T x(t) = 0, so a£(r) = 0 for all t. Such constant 
offsets have no effect on the Allan variance, which is 
non-zero only if the frequency difference of the two 
oscillators is changing with time. 

(ii) L inear frequency drift : 

x(r) = jet 2 (so y(t) = ct) 

In this case A 2 T x(t) = ct 2 , and 


1 Also, D. Percival, letter to J. Barnes, Aug. 31, 1982. 
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A linear frequency drift causes the Allan deviation 
to be proportional to r. Section IV shows one way 
of estimating the drift and removing it from the 
measurements. 


looks like a staircase, inspect your raw x or y data for 
periodic contamination. 

C. The Classical Random Clock Noise Model 


(iii) Higher powers of t . If n > 3 andx(r) = cf 1 9 then 

A 2 x(r) = cn(n - l)r 2 f" -2 + terms of lower degree 


Since this grows like a positive power of t , the Allan 
variance does not exist. (The limit in (3) is +«>.) 

(iv) A single frequency spike. This is the same as a step in 
phase or time. Suppose that x(t) jumps by an amount 
X Q during a measurement of duration T . Given r = T/m 
there is an index k such that A T x(kr) = X 0 and the 
other A t x(jt) are zero. Then 


A 2 r x(kt)=X Q , A*x((k + 1)/) = -X Q 


and the rest are zero. It follows from (4) that 


Sy(r, m) 


r(r-r) 


Thus, for r « T, the estimated Allan variance is 
approximately proportional to 1/r. A single frequency 
spike mimics white frequency noise (see below). This 
is appropriate, since a Poisson train of spikes is a form 
of white noise. 

(v) Periodic disturbances . These can be caused by daily 
temperature variations, periodic weather fluctuations, 
problems with the measurement system, or even prob- 
lems with the oscillators themselves. Their effect on 
Allan variance measurements can be bizarre. For exam- 
ple, let 


This is a combination of five “power-law” random proc- 
esses, each of which contributes its signature to the Allan 
variance. The model is specified by the one-sided spectral den- 
sity S y (f) of the fractional frequency j(t): 

s y (f) = h 2 f 2 + h/+ h 0 + hj - 1 + h_ 2 r 2 

- t V* 

a = - 2 

Table 2 gives the Allan variances of the five components of 
the model (Ref. 1). 

For example, if one sees o y {r) proportional to l/y/r over 
two decades or more, one usually interprets this as the effect 
of white frequency noise. You had better eyeball the x or y 
data versus time to make sure that the \\\fr is not caused 
by a single monster phase jump (see item (iv) above). The 
inference of a model from o y (f) is unreliable without some 
common-sense checks. 

For white phase and flicker phase, the o*(t) formulas are 
approximations that require 2 nf h T » 1. Moreover, the white 
phase formula can be extended to the case of a stationary 
time deviation process x(t) having an autocovariance function 
R x (r) = E x(t)x(t + r). Let = ^(0), the time variance. If 
there is a r 0 such that \R x (r)\ «R X ( 0) for t > r 0 , then 

3a 2 

°y^ 88 ~T > T>T o ( ? ) 

T 


x(r) = X Q cos (2 7 tvt + Q) 


III. Zero-Crossing Counter Technique 


Some computation yields 
o y ir) = 

j (y/2\cosQ\, if 2 vt = integer 

flZJsin 2 (t r*r)- j 

( 1 , otherwise 


This is the theoretical Allan deviation (passing to the 
limit in (4)). For 6 = 0, this function looks like Fig. 3. 
If your Allan deviation plot has a lot of wiggles or 


Figure 4 shows how the frequency and timing groups in 
the Communications Systems Research Section measure the 
time deviation x(r) of two oscillators. The oscillators are off- 
set in frequency by v b , which, in the present setup, is at most 
1 Hz. Oscillator 2 runs at the higher frequency. If oscillators 1 
and 2 have phases and 0 2 , then the beat-note signal is 

sin (27rjy-0(f)) 

where = 0j(O - 0 2 (O- zero-crossing detector triggers 
the counter at the beat-note upcrossing times t 0 , f 2 , *••. 

The counter, running continuously at frequency v c (1 MHz in 
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the present setup), records discrete approximations ?q, t [ , 
t ' 2 , ••• for further digital processing. The quantization errors 
tj- tj will be dealt with later. 

Let r b = 1/^, the nominal beat period. Then tj - ~ r 6 . 

Note that r b is not precisely defined by the data. It could be 
t Y - t 0 , or the average of all the tj - f ; -_ t over the duration of 
the test. Changing r b also changes 0(0, but only by const • f, 
so that frequency stability is not changed. 

The time deviation of the beat note is 


ity of the beat note must be no worse than 1%. The beat note 
itself must be a reasonably good clock, but not so good that its 
fluctuations are hidden by counter quantization. In most of 
our tests, this stability is 10~ 4 or better. There was one test, 
however, in which the beat periods tj - fluctuated by a 
factor of two or more. In this situation, the Allan variance esti- 
mate given below is invalid. It is thought, however, that there 
might be a more sophisticated algorithm that would still be 
able to estimate Allan variance for r » T b . 

Since the second differences of / are zero, we have 


wa - m - urn - 

x b Kl) 2w„ v. 2 nv. 

D D V 



v 0 A r V 2t U-i)r + t U-2)^ 


where x(f) is the time deviation of the oscillator pair. (Note 
the sign reversal.) Since the factor vjv b is large, typically 
10 6 to 10 8 , the fluctuations of the oscillators are magnified 
so that they can be measured. The Allan deviation of the oscil- 
lators is v b /v o times the Allan deviation of the beat note. 

The computation of the Allan variance of the beat note is 
not totally straightforward. By definition, the / th upcrossing 
tj occurs when the total phase of the beat note is 2ir]\ i.e., 

2? rjy y -0(L) = 271/ 

At time tj the time deviation of the beat note is 


Using this in (4) in place of A \x(jr) gives the Allan variance 
estimate 


S 2 
y 




where r = rr 0 and the duration of the test is mr. Look again 
at the numerical example in Section II. The x(f) column of 
Table 1 actually gives tj r in microseconds (r = 256,/ = 0 to 8, 
r b = 1 s) with the gross linear part jrr b already subtracted off. 
Since v b = 1 Hz, v Q = 10 8 Hz, the scale factor is 10~ 6 v b lv 0 

= 10- 14 . 




2irv v 


= jr b - t p /' = 0, 1,2, 


(another sign reversal), and hence the time deviation of the 
oscillators is 


Quantization Error. Let the counter frequency be v c . (Ours 
is 1 MHz.) My model for the recorded time tj of the / th up- 
crossing tj is 


V = t i +q i 


Now suppose that the aim is to measure o y (T) for r = rr b , 
r a positive integer. This could be done if we had the sequence 

*('<))’ x ( r o + T )> x ( { 0 + 2t ^> “ ■ 

What we do have is the sequence 

^ 0 ). X (t r l X ( ( 2 ( 9 ) 

which is x at times that are not exactly spaced by r. This is all 
right as long as the periods t /r - differ from r by at 

most 1%, say. In other words, the fractional frequency stabil- 


where tj is a multiple of l/i> c , and q f , the quantization error, 
satisfies 0 < < 1 jv c . As an approximation, assume that the 

processes tj and Qj are uncorrelated, and that is a process of 
independent, uniform random variables. This approximation 
also yields Sheppard’s second-moment correction in statistics. 
Then the alone contribute an Allan variance 



(This follows from (7).) For example, if = 5 MHz, v b ~ 1 Hz, 
v c = 1 MHz, then a q (j) = 10” 13 /r. One should a^(r) 

from the measured a£(r) (from tj) to get the true aj(r) (from 
tj). In other words, a^(r) contributes to the measurement- 
system noise floor. 
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IV. Drift Removal 

During a long test of a pair of hydrogen masers, it is often 
apparent that the frequency deviation y(t) appears to be dom- 
inated by a linear component, which also causes o y (r) to 
increase like r for large r, thus masking the effect of the ran- 
dom fluctuations. See Fig. 5 and the upper curve of Fig. 6. 
The model is 

x(0 = a 0 +a x t + ^ct 2 +x 0 (t) 


In other words, drift rate equals average frequency at the end 
of the run, minus average frequency at the beginning of the 
run, divided by the time span between the midpoints of the 
averaging intervals [0, r c ] and [T - r c X\ * Take 

T 

T c 6.29 

a value that minimizes the variance of c 2 in the presence of 
flicker-frequency noise. This value is not critical. In terms of 

m. 


or 

<2*o 

*') = a i +Ct + ^F 

where c is the frequency drift rate (A/// per second) and 
x Q (t) is a mean-zero process. The constants tf 0 , a x are irrele- 
vant for a frequency stability study. Let a 2 (r), o 2 0 (r) be the 
Allan variances of x 0 (f). We would like to estimate c and 
to use the estimate c to extract an estimate of aj 0 (r), the 
Allan variance of the underlying fluctuations. 

A, Estimating 8 From x(f), 0 < t < T 

Assume a sample time of r 0 . The c estimators given below 
are all unbiased. 

Method 0. Least-squares quadratic fit to x(r): 

*(0 = a Q + a 1 t+~c Q t 2 

where a 0 , a j , c 0 are chosen to minimize 2 [£(7) - x(r)] 2 . This 
method is optimal (minimal variance) in the presence of white 
phase noise. 


x(7) - X (T- t ) - x(t c ) + jc(0) 

= (10) 

C K c ' 

so that only four values of x(f) are needed. 

This estimator, although suboptimal for all the classical 
noises, performs well in the presence of white, flicker, and ran- 
dom walk frequency noises. If * 0 (0 is random walk fre- 
quency, then the standard deviation of c 2 is 

a vo( T C ) 

a(c 2 ) = 4.6-^p- (11) 

and this formula is pessimistic for white and flicker frequency. 
Note that you have to remove drift in order to estimate 
v y o(T c ); see (13) below. 

B. Removing Drift 

Suppose that some unbiased estimate c is obtained, perhaps 
one of the foregoing, perhaps some other. An estimate of the 
residual x 0 (r) is 

x 0 (t) = x(t)-jct(t-T) 


Method 1. Least-squares linear fit tojv(f,T 0 ): 

9(0 = «!+?! t 

where , c x are chosen to minimize 2 [.v’(f') -7(f,T 0 )] 2 ■ This 
method is optimal in the presence of white frequency noise. 

Method 2. Mixed second difference of phase. This is the one 
that has been used here for tests of frequency standards (Refs. 
4 and 7). No one else is thought to have used it. So far, no 
complaints have been received. Let 

„ 9(T,r c )-y(T c ,T c ) 

C 2 T-T 

C 


See Fig. 7 for the effect of this (c = c 2 ) on the data of Fig. 5. 
Then 

A** 0 (0 = A 2 T x(t)-ct 2 
is used in the estimate 

m r -1-7 

s >- m) 

for a 2 0 (r). Expanding the square gives 

S 2 y0 (r,m) = S 2 y (T,m)- T — c 2 t +^-(c t -c) 2 (12) 
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where 


x(T *) - x(T - r) - x(t) + Jt(0) 


and T- mr. If c = c 2 , then 


(2) Relative linear frequency drift is negligible or removed. 
In other words, the second differences of x A , x B , x c 
have mean zero. 

If the assumptions are satisfied, then, for a fixed r, the 
processes A 2 x a , A 2 x B , A 2 x c are orthogonal. Therefore, 


O 3 ) 

is used for estimating o 2 q (t c ), needed in (11). The lower 
curve in Fig. 6 is S y0 (T,m ), the upper curve is S y (r,m), and 
the straight line is (6) with c = c 2 . 

CAUTION 

If c is based on the test data x(t) (0 < t < T), T = mr, 
m is small (<5, say), and random walk frequency noise 
dominates x 0 (t), then S 2 0 (r y m) severely underestimates 
a vo( T ) on a ve rage, and has a larger variance than 
Sy(T,m). If c = C 2 , the bias of S 2 0 is given by Ref. 7. 

m 10 7 5 4 3 2 

Bias -16% -24% -34% -43% -59% -89% 

The error bars in Fig. 6 are the naive ones, and do not 
reflect these biases. The lesson here is to use a long-term 
estimate of c if possible. For example, cavity retuning of 
a hydrogen maser gives an independent measurement of 
frequency drift. 


V. Three-Cornered Hat 

The previous material has been about the relative stability 
of a pair of oscillators. If one has a triplet of oscillators, one 
may be able to estimate the stability of each one individu- 
ally by using the three pair-comparisons. The setup, as shown 
in Fig. 8, has oscillators A, B, C, and pair-channels 1, 2, 3. 
We measure the pair Allan variances a](r), o\{r), o 2 (t) and 
would like to compute o 2 (t ), a|(r), a£(r) for the individual 
oscillators. 

Assumptions 

(1) The phase fluctuations of the three oscillators are 
independent. 


and so 


o] = °l + a l 


°\ = n l + °A 


°l = ° 2 A +a l 


° 2 a 

°l = ^ (°l + ~ °l) 
a l = i k + ■ °l) 


(14) 


(15) 


Figure 8 shows the geometric interpretation. 

Difficulties 

(1) Since frequency standards are sensitive to changes in 
the environment (temperature, pressure, humidity, 
magnetic field), it seems that frequency standards in 
the same room must violate assumption 1 for long-term 
fluctuations, say, for 1000 s or more. 

(2) The pair-channel o 2 (t) values are not available, just 
estimates such as iS 2 (r,m). When estimates are substi- 
tuted for the a 2 values in the right side of (15), it often 
happens that some oscillator o 2 comes out negative. 
One interpretation of this is that a confidence interval 
for that a goes all the way down to zero. For example, 
if *S^ > S\ + S^, then the o 2 A estimate is negative. (If 
a triangle has sides S 19 S 2 , S 3 , then the angle between 
S 2 and S 3 is obtuse.) If the stronger inequality S x >S 2 
+ S 3 holds, then something is wrong with the measure- 
ments. (There is no triangle with sides S x , S 2 ,S 3 .) 
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Table 1. Data from a stability test of two hydrogen masers 


t 

x(t) 


AjX(r) 

& 2 x(t) 


v (,) 

A 3*('> 

0 

0 







1 

658 

658 






2 

1229 

571 

-87 

1229 




3 

1701 

472 

-99 



1701 


4 

2333 

632 

160 

1104 

-125 

1675 


5 

2991 

658 

26 



1762 


6 

3493 

502 

-156 

1160 

56 

1792 

91 

7 

4095 

602 

100 



1762 

87 

8 

4690 

595 

-7 

1197 

37 

1699 

-63 


Table 2. Spectral density and Allan variance of the components 
of the classical clock noise model 


a 

Name 

V n 

%( t ) 

2 

White phase 

h 2 f2 ’ f< f h 

h 1*h .22 
4 7T T 

1 

Flicker phase 

V’/<4 

h, -4-r In (8.88 f r) 
1 4tt r 

0 

White frequency 

V f h = + ~ 


-1 

Flicker frequency 


/i In 4 

-2 

Random walk 
frequency 


, 27r 2 r 

h -2 3 




RUN: 99-1 A 
START: 1230,12-23-80 
END: 0820, 12-31-80 
CHAN 3: DSN2-DSN3 
FCEN: 1.000+08 Hz 
REMOVED DRIFT: 0.000/day 
SAMPLE TIME: 4.094+03 s 




Fig. 2. Estimated Allan deviation of a pair of hydrogen masers 
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Fig. 3. Allan deviation for x(f) = cos (2irvt) 


OSCILLATOR 1 


OSCILLATOR 2 


Fig. 4. Frequency stability test setup 
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TIME, 10° s 

Fig. 5. x and / before drift removal 


RUN: 99-1 A 

COMPUTED DRIFT: 


START: 1230, 12-23-80 

-5.179-14/day 


END. 0820,12-31-80 
CHAN 1: DSN2-DSN1 

+6.678-15 


FREQ: 1 .000+08 Hz 
BW: 1.000+00 Hz 
TSEP: 7.993+00 s 

BEFORE — v 



^ \y 



y 



81 

s 

V 40 


' — AFTER 



Fig. 6. Estimated Allan deviation before and after drift removal 











